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We study the survival of charmonium states in a strongly-coupled quark-gluon plasma in the 
presence of dissipative forces. We consider first-order dissipative corrections to the plasma equation 
of motion in the Bjorken boost-invariant expansion with a strongly-coupled equation of state for 
QGP and study the survival of these states with the dissociation temperatures obtained by correcting 
the full Cornell potential not its Coulomb part alone with a dielectric function encoding the effects of 
deconfined medium. We further explore the sensitivity of prompt and sequential suppression of these 
states to the shear viscosity-to-entropy density ratio, rj/s from perturbative QCD and AdS/CFT 
predictions. Our results show excellent agreement with the recent experimental results at RlllC. 
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I. INTRODUCTION 

One of the amazing discoveries of experimental 
measurements at RHIC is the surprising amount of 
collective flow exhibited by the outgoing hadrons 
which is evidenced in both the single-particle trans- 
verse momentum distribution [l| and the asymmet- 
ric azimuthal distribution around the beam axis , 
quantified by the parameters, vi and V2, respectively. 
The parameter, V2, in particular, was expected to be 
much smaller at RHIC than it was at the lower ener- 
gies of the SPS 3 ; in fact, it is about twice as large. 
Theoretical calculations cannot predict such higher 
values of vi and V2 unless the partonic cross-sections 
are artificially enhanced by more than an order of 
magnitude over perturbative QCD predictions [J. 
Thus quark-gluon matter created in these collisions 
is strongly interacting (SIQGP), unlike the type of 
weakly interacting quark-gluon plasma (wQGP) ex- 
pected to occur at very high temperatures on the 
basis of asymptotic freedom [^. On the other hand, 
lattice QCD calculations yield an equation of state 
that differs from an ideal gas only by about 10% 
once the temperature exceeds l.STc @. Further- 
more, perfect fluid dynamics with zero shear and 
bulk viscosities reproduces V2 well up to transverse 
momenta of order 1.5 GeV/c; however, it cannot be 
zero. Indeed, the calculations within AdS/CFT sug- 
gests that T] > s/{4tt). 

Charmonium (J/V") suppression has long been 
proposed as a signature of QGP formation^] 
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and has indeed been seen at SPS [I and RHIC 
experiments[l|. The heavy quark pair leading to the 
J/ip mesons are produced in nucleus- nucleus colli- 
sions on a very short time-scale ^ l/2mc, where 
rric is the mass of the charm quark. The pair de- 
velops into the physical resonance over a formation 
time and traverses the plasma and (later) the 
hadronic matter before leaving the interacting sys- 
tem to decay (into a dilepton) to be detected. This 
long 'trek' inside the interacting system is fairly 'haz- 
ardous' for the J/tjj- Even before the resonance is 
formed it may be absorbed by the nucleons stream- 
ing past it 0]. By the time the resonance is formed, 
the screening of the colour forces in the plasma may 
be sufficient to inhibit a binding of the cc . The 
resonance(s) could also be dissociated either by an 
energetic gluon [lOl . ITll | or by a comoving hadron. 
A complete and deep understanding of the various 
stages of absorption is still missing and is under 
intense investigations [T^ . In order to disentangle 
these different effects, we must know the properties 
of quarkonium in medium and determine their dis- 
sociation temperatures. 

The physics of a given quarkonium state is en- 
coded in its spectral function [Tsl. [T^ . So following 
how the spectral function changes with temperature 
can give us a theoretical insight to the temperature- 
dependence of quarkonium properties. There are 
mainly two lines of theoretical studies to determine 
quarkonium spectral functions: the first one is the 
potential models [H, O, [13, [Tl] which have been 
widely used to study quarkonium, but their applica- 
bility at finite temperature is still under scrutiny and 
the second one is the lattice QCD[1^ which provides 
the most straightforward way to determine spectral 
functions, but the results suffer from discretization 
effects and statistical errors, and thus are still incon- 
clusive. These two approaches show poor matching 
between their predictions because of the uncertain- 
ties coming from a variety of sources. None of the 
approaches give a complete framework to study the 
properties of quarkonia states at finite temperature. 
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However, some degree of qualitative agreement had 
been achieved for the S-wave correlators. The find- 
ing was somehow ambiguous for the P-wave corre- 
lators and the temperature dependence of the po- 
tential model was even qualitatively different from 
the lattice one. Refinement in the computations of 
the spectral functions have recently been done by 
incorporatin g th e zero modes both in the S- and P- 
channels [20l |21|. It was shown that, these contribu- 
tions cure most of the previously observed discrep- 
ancies with lattice calculations. This supports the 
use of potential models at finite temperature as an 
important tool to complement lattice studies. 

It is now pertinent to ask how a quarkonium be- 
haves in strongly interacting QGP unlike in wQGP. 
The propagation of charm quark in SIQGP is dif- 
ferent from wQGP because in SIQGP there will be 
a rapid equilibration due to the multiple momen- 
tum exchange in the momentum space but in the 
coordinate space, the motion is diffusive in nature 
and slower |23|. Therefore, the usual picture of char- 
moniun suppression Q may not be true in SIQGP. 
Moreover, in the RHIC era (small /is), it has been 
confirmed from the recent studies that the transi- 
tion from nuclear matter to QGP is not a phase 
transition, rather a crossover |23l|. It is then rea- 
sonable to assume that the string-tension does not 
vanish abruptly at the deconfinement point, so we 
may expect presence of non-pcrturbative effects such 
as non-zero string tension in the deconfined phase 
and one should study its effects on heavy quark 
potential even above Tc- This is indeed compat- 
ible with the strongly interacting nature of QGP 
where large number of non-pcrturbative effects are 
expected. This issue, usually overlooked in the liter- 
ature where only a screened Coulomb potential was 
assumed above Tc and the linear/string term was 
assumed zero, was certainly worth of investigation. 
This was exactly we had done recently [23| where 
a medium-modified heavy quark potential was de- 
rived by correcting the full Cornell potential not its 
Coulomb part alone with a dielectric function encod- 
ing the effects of the deconfined medium. We found 
that this approach led to a long-range Coulomb po- 
tential with an (reduced) effective charge in addi- 
tion to the usual Debye-screened form. With such 
an effective potential, we investigated the effects of 
pcrturbative and non-pcrturbative contributions to 
the Debye masses on the dissociation of quarkonium 
states. Finally we determined the binding energies 
and dissociation temperatures of the ground and the 
first excited states of charmonium and bottomonium 
spectra |2l,[2|. 

Now let us consider a central collision in a nu- 
cleus nucleus collision, which results in the forma- 
tion of quark gluon plasma at initial time . Let us 
concentrate on the region of energy density greater 
than a characteristic energy density, say, screening 
energy density Cs which encloses the plasma which 
is dense enough to cause the melting of a particu- 



lar state of quarkonium. We assume the plasma to 
cool, according to Bjorken's boost invariant (longi- 
tudinal) hydrodynamics. If the plasma has cooled to 
a energy density less than es, the cc pair would es- 
cape and quarkonium would be formed. If however, 
the energy density is still larger than e^, the reso- 
nance will not form and we shall have a quarkonium 
suppression. It is easy to see that the depen- 
dence of the survival probability will depend on how 
rapidly the plasma cools. If the initial energy den- 
sity is sufficiently high, the plasma will take longer 
to cool and only the pairs with very high will 
escape. If however the plasma cools rapidly, then 
even pairs with modera te will escape. There are 
many attempts [H, [21, [2^ where the above picture 
of plasma expansion dynamics had been employed to 
study charmonium suppression, but in their works 
following points arc missing: a) the effects of dissi- 
pative forces was not included in hydrodynamic ex- 
pansion. The presence of dissipative term in the 
energy-momentum tensor [30l | causes the plasma to 
expand slowly and the system will take longer to 
reach (screening time r^) and only the pairs with 
very high will escape resulting overall more sup- 
pression, b) The EoS employed in their works was ei- 
ther ideal or bag model EoS which was then directly 
used to calculate two crucial factors in quantifying 
the suppression in relativistic heavy ion collisions: 
the energy density corresponding to the dissocia- 
tion temperature (To) and the time elapsed by the 
system (ts) to reach Cg by the expansion (cooling) 
of the system. As we now know the matter formed 
at RHIC is far from its ideal limit even at T > Tc so 
the EoS employed in their works was not reliable to 
capture the properties of strongly-interacting mat- 
ter. This is evident from the substantial difference in 
the numerical values of the crucial quantities calcu- 
lated in the ideal and strongly-interacting equation 
of state (in the Tables II- V). 

The main motivation of this article is to rem- 
edy the above shortcomings in multifold respects: i) 
First we use an appropriate equation of state (EoS) 
which should reproduce the lattice results verifying 
the strongly-interacting nature of QGP. This also 
explains the non-ideal behaviour of QGP through 
the speed of sound which will play an important 
role in the expansion dynamics, ii) Then we study 
hydrodynamic boost-invariant Bjorken expansion in 
(1 + 1) dimension with the EoS discussed in i) as an 
input. In addition we explore the effects of dissipa- 
tive terms on the hydrodynamic expansion by con- 
sidering the shear viscosity 77 up to first-order in the 
stress-tensor. Basically we consider the ratio of the 
shear viscosity-to-entropy density (jy/s) as 1/A-k and 
0.3 which was predicted from the AdS/CFT corre- 
spondence [U and pcrturbative QCD [32| , respec- 
tively. We also consider ideal hydrodynamics with 
77/s = for the sake of comparison, iii) The most 
important point is to know the properties of quarko- 
nium in the medium. For that we need to know the 
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correct form of medium-modified potential to deter- 
mine the dissociation temperatures of charmonium 
states in a static, thermal medium. For this pur- 
pose we follow our recent work where medium mod- 
ified potential [23, was obtained by correcting 
the full form of the Cornell potential through a di- 
electric function embodying the medium effects, iv) 
Finally we study the survival of charmonium states 
with all the refinements discussed above. The ex- 
tent of suppression will be decided by a competition 
among the dissociation scale, the transverse momen- 
tum of J/ip^s and the rate of expansion, making it 
sensitive to such details as the speed of sound j 26, 2^ 
through the EoS. Thus a study of the survival of J/ tp 
is poised to provide a wealth of information about 
the dissociation mechanism of cc states in a hot QCD 
medium, the evolution of the plasma and its prop- 
erties. 

The paper is organized as follows. In Sec. II. A, 
we briefly describe the strongly-interacting equation 
of state developed by Bannur [Sl, [33| and deter- 
mine the pressure, energy density and the speed of 
sound etc. In Sec. II. B, we then use the aforesaid 
equation of state to study boost-invariant (1+1) di- 
mensional longitudinal expansion in the presence of 
viscous forces. In Sec. III. A, we study the dissocia- 
tion phenomenon of J/ ijj by studying the in- medium 
modifications to the heavy quark potential and its 
suppression in a longitudinally expanding QGP in 
Sec. III. B. Results are presented in Sec. IV with a dis- 
cussion on important aspects of our approach. Fi- 
nally in Scc.V, we present the conclusions and future 
prospects of the present work. 



II. STRONGLY-INTERACTING QGP AND 
LONGITUDINAL EXPANSION OF QGP 

In this section, we first discuss the equation of 
state for SIQGP. We shall then discuss hydrodynam- 
ical expansion in the presence of viscous forces with 
the SIQGP EoS as an input. 



A. Equation of state 

The equation of state for the quark matter pro- 
duced in RHIC is a very important observable and 
the properties of the matter are sensitive to it. The 
expansion of QGP is quite sensitive to the equation 
of state(EoS) through the speed of sound. Both ex- 
perimental |35[ and theoretical (lattice) results [36j 
show that matter formed near or above Tc is non- 
ideal. There have been many attempts to explain 
such a matter using various models such as bag mod- 
els, other confinement models [33], quasi particle 
models [3711 . strongly interacting quark gluon plasma 
(sQGP) [23 etc. 

Here we propose that the QGP near Tc is in fact 
what is called strongly coupled plasma [S^l, widely 



studied in QED plasma where the plasma parame- 
ter, F, defined as the ratio of the average potential 
energy to the average kinetic energy of the particles, 
is of the order of 1 or larger. An extensive study 
of strong-coupled plasma in QED with proper mod- 
ifications to include colour degrees of freedom and 
the strong running coupling constant gives an ex- 
pression for the energy density as a function of the 
plasma parameter, F 



(2.7 + Uex(F))nr 



(1) 



where the first term corresponds to the ideal EoS 
and the second term, Uox(F), gives the non-ideal (or 
excess) contribution to EoS as 



(2) 



So, the scaled energy density (in terms of Stefan- 
Boltzmann limit) is given by 



2.7 2 



e(r) = — = 1 - — ^F^/^ , (3) 

where esB = 3a fT'^ with degrees of freedom aj 
(16 + 21ny/2)7r^/90 and using the relation, e 

dp 
dT 



— we get the pressure 



P_ 



dTT'^e{V{T)) /T^ 



(4) 

where po is the pressure at some reference temper- 
ature Tq. If we calculate p{T)/T'^ versus T for 
pure gauge, 2-flavor and 3-flavor QGP, a surprisingly 
good fit with the lattice data was found [33|, [131 • 

Another important observable, the speed of sound 
which appears in the equation of motion for the ex- 
pansion through the relation. 



dp 
'd'e' 



(5) 



It is found that strongly-interacting matter created 
at RHIC has a signature of non-zero interaction mea- 
sure, A = e — 3p which measures the deviation of the 
EoS from its ideal counterpart. This quantity has 
been extensively studied in lattice QCD (3^ which 
shows that A > and asymptotically approaches to 
the ideal value (zero). This implies that < 1/3 at 
T = Tc or even much above Tc- This is exactly found 
in Bannur's model [33| where is very much less 
than its ideal limit around and increases towards 
its ideal value. In view of the excellent agreement 
with lattice simulations the above phcnomcnologi- 
cal EoS is a right choice for the strongly-interacting 
matter possibly formed at RHIC to calculate the 
thermodynamical quantities viz. screening energy 
density [eg), the speed of sound etc. and also to 
study the hydrodynamical expansion of plasma. 
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B. Longitudinal expansion in the presence of 
dissipative forces 

The energy momentum tensor of the plasma in the 
absence of dissipative forces is written as: 



T^"' = (e + p)u^'u'' + g^'•'p 



(6) 



If the effects of viscosity are neglected, conservation 
of energy-momentum tensor leads to 



(7) 



We now consider the Bjorken boost invariant longi- 
tudinal flow in (1+1) dimension where the equation 
([7|) simplifies to 



dre = — 



e+p 



(8) 



The effect of the speed of sound is seen immedi- 
ately through the EoS: p = c^e as 



e(r)T 



const. 



(9) 



With the condition e ^ T'^ and taking the speed of 
sound = 1/3, the above equation translates to the 
cooling law: 



(10) 



where Ti is the initial temperature and e(Ti) is the 
initial energy density. 

Now, we study the correction to the equation of 
motion ([7]) in the presence of dissipative forces in the 
form of shear viscosity, rj which is itself a function 
of temperature. The hydrodynamic evolution pre- 
serves the boost-invariance, which in ultrarelativistic 
heavy ion collisions is expected to be realized near 
mid-rapidity [40l |. In the presence of viscous forces 
the energy- momentum tensor is written as. 



T^"' = (e + p)u''u'' + g'^'p + ir"" , 



(11) 



where tt^'' is the stress-energy tensor. Dissipative 
corrections to the first-order in the gradient expan- 
sion are recovered by setting the relaxation time to 
zero. This leads to: 



(12) 



where 77 is the co-efficient of the shear viscosity and 



is the symmetrized velocity gradient. 



In (1-|-1) dimensional Bjorken expansion in the 
first-order dissipative hydrodynamics, only one com- 
ponent tt'^'^ of the viscous stress tensor is non-zero. 
In this case the equation of motion reads. 



Ore = h 7-^ 



477 
3^2 



(13) 



The solution of the above equation in the case of 
constant value of rj/s is known analytically [4ll 
and is given by, 



T(t) - 



1/3 



1 



27y 



SsTiTi 



1 



2/3 



(14) 

The first term in the RHS is the same as in the 
case of zeroth-order (non-viscous) hydrodynamics 
and the second term is the correction arising from 
constant rj/s which causes the system to expand 
slowly compared to the perfect fluid 77 = pO|) . It 
is important to note that temperature will be max- 
imum after a span of time 



- 1 S TjTi 

'^111 ax I o n 

V 2 



-3/2 



(15) 



For times t > T,„ax, the temperature decreases 
with time, as expected for matter undergoing ex- 
pansion. For early times r < Tmax, however, the 
solution shows an unphysical reheating which 
will not be discussed here. 

We shall now study the solution of Eq. (fT5|) for 
the strongly-interacting equation of state discussed 
in Sec. II. A where is always less than 1/3 and is a 
function of temperature. The solution of Eg. p^ is 
obtained as, 



3t^ 



on 



const , 



(16) 



where a = (^) Tf r, 
(1 - cI)t^ and (1 - 



and 



and f~ are given by 

respectively. The first 
term in both LHS and RHS accounts for the con- 
tributions coming from the zeroth-order expansion 
and the second term is the first-order viscous cor- 
rections (In other words, Eq. (fT6|) reduces to Eq.®) 
for T] = 0.) 

In our work wc consider three values of the shear 
viscosity-to-entropy density ratio to see the effects of 
nonzero values of the shear viscosity on the expan- 
sion. The first one is from perturbative QCD [3^ 
calculations where 77/s is =0.3 near T ~ Tc — 2Tc. 
The second one is from AdS/CFT studies [3l| where 
77/s = l/47r (^ 0.08). Finally we consider 77/5=0 (for 
the ideal fluid) for the sake of comparison. Wc shall 
employ Eq. (jl6[) to study the charmonium suppres- 
sion in a strongly interacting QCD medium formed 
in a ultra-relativistic heavy-ion collisions in the next 
section. 



III. SUPPRESSION OF J/tl^ IN A 
LONGITUDINALLY EXPANDING PLASMA 

Let us now turn our attention to the most impor- 
tant quantity, viz., the J/^p survival probability. Re- 
call that, to study the fate of charmonium in SIQGP, 
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we need to consider (i) the equation of state for 
SIQGP, (ii) the effects of dissipative forces (shear 
viscosity) on the expansion, and (iii) an appropri- 
ate criterion for the dissociation of charmoniun in 
hot QCD medium. The first two issues have ah'cady 
been discussed in the previous section. Here, we first 
discuss the dissociation of charmonium states by giv- 
ing a brief derivation of medium-modified potential 
in Sec. III. A. We also discuss the nature of dissoci- 
ation by examining the effects of perturbative and 
non-perturbative terms in the Debyc masses. Finally 
we proceed to derive the survival probability of J /ip 
in an expanding SIQGP in Sec. III. B. 



A. Dissociation of charmonium in SIQGP 

The large distance property of the heavy quark 
interaction is important for our understanding of 
the bulk properties of the QCD plasma phase, 6.17. 
the screening property of the quar k g luon plasma 
[23, H^l , the equation of state |46|. 1471' and the or- 
der parameter (Polyakov loop) [48l. I49I. ISOl [5]| . In 
all of these studies deviations from perturbative cal- 
culations and the ideal gas behavior are expected 
and were indeed found at temperatures which are 
moderately larger than the deconfinement temper- 
ature. This calls for quantitative non-pcrturbative 
calculations. Recent findings [2^ indicate the phase 
transition in full QCD appears as an crossover rather 
than a 'true' phase transition with the related sin- 
gularities in thermodynamic observables. In light of 
this finding, one can not simply ignore the effects 
of string tension between the quark-antiquark pairs 
beyond Tc- This is indeed a very important effect 
which needs to be incorporated while setting up the 
criterion for the dissociation of quarkonia in QGP. 

In our recent works [13, [2^, this issue has success- 
fully been addressed for the dissociation of quarko- 
nium in QGP. In these works, an appropriate form 
of the medium modified inter-quark potential have 
been obtained. The authors have considered the 
Cornell form of the potential: V{r) = —a/r + ar 
and correct its Fourier transform (FT) V{k) as 



V(k) 



Vjk) 
e{k) 



(17) 



where e{k) is the dielectric permittivity given in 
terms of the static limit of the longitudinal part of 
the gluon self-energy [52 : 



eik) 



1 



nL(o,A:,r) 



1 



k^ 



(18) 



and V{k) is the Fourier transform (FT) of the Cor- 
nell potential: 



V(fc) = -v/(2M 



4cr 



k^ y^nk^' 



(19) 



Substituting Eqs.([TH]) and ([HD into ([T7D and 
then evaluating its inverse FT one obtains the r- 
dep endence of the medium modified potential [2^ . 



-mDr) 



mj~) r 

2a 2a 

— — 9 \ aniD , 



(20) 



where the constant terms are introduced to yield 
the correct limit of V{r, T) as T ^ 0. Such terms 
could arise naturally from the basic computations 
of real time static potential in hot QCD|5^ and 
from the real and imaginary time correlators in a 
thermal QCD medium [m|. The medium modified 
form of the potential thus obtained has an additional 
long range Coulomb term with an (reduced) effec- 
tive charge in addition to the conventional Yukawa 
term. A new picture of quarkonium suppression has 
emerged out of these studies where both the range 
and the charge of the potential get screened due 
to in-medium modifications. The binding energies 
and dissociation temperatures for various quarko- 
nium states have been determined by solving the 
Schrodingcr equation numerically with the poten- 
tial (|20|) . We considered three possible forms of the 
Debye masses (mj^) in the potential: the leading- 
order term in QCD coupling [in^^ ) which is a per- 
turbative result [2^ . the non-perturbative correc- 
tions to the leading-order {m^) [5^ and the lattice 

parametrized form {mj^ = 1.4m^^) [56l | to exam- 
ine the effects of perturbative and nonperturbative 
contributions on the dissociation mechanism. The 
upper limit of the dissociation temperature of a par- 
ticular state is obtained when its binding energy is 
of the order of temperature of the system: i^bin ~ T 
so that the state becomes feebly bound and thermal 
fluctuations can excite it to the continuum. The 
perturbative result for the Debye mass in the po- 
tential gives rise higher values (in Table II) of the 
dissociation temperatures compared to the spectral 
function studies [H, [13] calculated in a potential 
model but smaller than the results from the anal- 
ysis of lattice (temporal) correlators [11] (in Table 
IV). On the other hand, when we use the lattice 
parametrized form (by fitting the lattice free-energy 
for color-singlet states) for the Debye mass, the dis- 
sociation temperatures (given in Table III) exactly 
matches (discussed in in details in Ref. [l^) with the 
results of the potential model calculations [5p, [53| 
but much smaller than the lattice studies [58|| (Ta- 
ble IV). Finally the non-perturbative corrections of 
0{g^T) and 0{g^T) to the perturbative result [11] in 
Debye mass give unrealistically smaller values which 
will not be considered in the present work. 

In addition, we also take the advantage to demon- 
strate the flavor dependence of the dissociation pro- 
cess. Hence we determined the dissociation tempera- 
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tures for the 2-flavor QGP (Table V) with the Debye 
mass only in the leading-order term in strong cou- 
pling. The dissociation temperatures are found to 
be higher than the 3-flavor case (Table II). However 
the detailed study of the flavor dependence was given 
in our earlier works [23, [2^ ■ In the present analy- 
sis, we have used the four sets of the dissociation 
temperatures (Tables II- V) to study the centrality 
dependence of survival probability of charmonium 
as a function of number of participants Np^^^ in an 
expanding (strongly interacting) QCD medium. 



B. Survival probability 

We now have all the ingredients to write down 
the survival probability and we closely follow Chu 
and Matsui for this. In the work of Chu and Mat- 
sui [2^, the px dependence of the survival proba- 
bility of charmonium was studied by choosing the 
speed of sound = 1/3 (ideal EoS) and the ex- 
treme value = 0. This work was further general- 
ized by invoking the various parameters for Au-Au 
collisions at RHIC in [23| to include the effects of re- 
alistic EoS in an adhoc manner by simply choosing a 
lower value of = 1/5. Instead of taking arbitrary 
values of wc tabulated the values of in Tables 
II-V corresponding to the dissociation temperatures 
calculated from Bannur model [33l | where the values 
of match perfectly with the lattice results. More- 
over, in the light of recent experimental finding from 
RHIC, one can not ignore the viscous effects while 
studying charmoniun suppression. Here, we address 
these issues. Note that in our analysis, we follow 
the idea of Chu and Matsui [2^ but with the more 
appropriate criteria for the charmonium dissociation 
and the dissipative hydrodynamic expansion of the 
plasma. 

Let us take a simple parametrization for the initial 
energy-density profile on any transverse plane : 







0{Ri 



-r) 
. (21) 

where r is the transverse co-ordinate and Rt is the 
transverse radius of the nucleus. One can define an 
average energy density < > as 



< Ci > = 



j 2TTrdre{r) 



so that 



£, = (1 +/?)(£,). 



(22) 



(23) 



We have taken f3 = 1, which will reflect the pro- 
portionality of the deposited energy to the nuclear 
thickness. The average energy- density is obtained 



TABLE L Kinematic characterization of Au-|-Au colli- 
sions at RHIC 01 



Nuclei ^/snn 


^part 


<e>i 


Rt 


(GeV) 




(GeV/fm^*) 


(fm) 




22.0 


5.86 


3.45 




30.2 


7.92 


3.61 




40.2 


10.14 


3.79 




52.5 


12.76 


3.96 




66.7 


15.69 


4.16 




83.3 


18.58 


4.37 


Au+Au 200 


103.0 


21.36 


4.61 




125.0 


24.38 


4.85 




151.0 


27.37 


5.12 




181.0 


30.52 


5.38 




215.0 


34.17 


5.64 




254.0 


37.39 


5.97 




300.0 


41.08 


6.31 




353.0 


45.09 


6.68 



from the Bjorken formula: 



1 



nRj^Ti 



dEr 
dy 



(24) 



where Et is the transverse energy deposited in the 
collision. The appropriate characterization of kine- 
matic quantities in Au+Au collisions is presented in 
Table I. 

The time r, when the energy density drops to the 
screening energy density is estimated from Eq. (|1( 
as 



^{r) 



4a 
3f? 



(25) 



where ei{r) = e{Ti,r) and is (1 — cDt^. The 
critical radius r^, is seen to mark the boundary of 
the region where the quarkonium formation is sup- 
pressed, can be obtained by equating the duration 
of screening Ts{r) to the formation time tp ~ JTp 
for the quarkonium in the plasma frame and is given 
by: 



RTil^A)^eil-A) , 



(26) 



where A is given by 

A = 



1 



1l 



4a 



1 4a 



(27) 



7 



with 



(1 



The quark-pair will escape 



the screening region (and form quarkonium) if its 
position r and transverse momentum are such 
that 



\r + TFPT/M\ > Ts. 



(28) 



Thus, if 4> is the angle between the vectors r and pr, 
then 

cos^ > [ir^,-r^)M-Tj,p'^j./M]/[2rTFPT], 

(29) 

which leads to a range of values of (j) when the 
quarkonium would escape. Now we can write for 
the survival probability of the quarkonium: 



r dr 



<^'max 



.P(r,pT) 



2tt 



r dr P{r,pT) 



I 

(30) 



where (/)max is the maximum positive angle (0 < 
< tt) allowed by Eq. 



(31) 



cos^^ \y\ if — 1 < y < 1 
if 2/ > 1 



where 



V 



'')M~Tlv\lM\l{2rTFVT\. (32) 



and P is the probability for the quark-pair produc- 
tion at (r, Pt), in a hard collision which may be 
factored out as 



P(r,pT) = /(r)5(pT), 
where we take the profile function f(r) as 

^2 



/(r) cx 



1 - 



B{Rt - r) 



(33) 



(34) 



with OL = 1/2. 

Often experimental measurement of survival prob- 
ability at a given number of participants ( A^p^^ ) '^i' 
rapidity (y) is reported in terms of the p^-integrated 
yield ratio (nuclear modification factor) over the 

range p^™' < Pt ^ P^^^ whose theoretical expres- 
sion would be 



{S{Pr)) = 



/min '^PtS{Pt) 

T 

pmax 



(35) 



The production of J/ip mesons in hadronic reac- 
tions occurs in-part through production of higher 



excited cc states and their decay into quarkonia 
ground state. Since the lifetime of different sub- 
threshold quarkonium states is much larger than 
the typical life-time of the medium which may be 
produced in nucleus-nucleus collisions so their de- 
cay occurs almost completely outside the produced 
medium. This means that the produced medium can 
be probed not only by the ground state quarkonium 
but also by different excited quarkonium states. 
Since, different quarkonium states have different 
sizes (binding energies), one expects that higher ex- 
cited states will dissolve at smaller temperature as 
compared to the smaller and more tightly bound 
ground states. These facts may lead to a sequential 
suppression pattern in Jf-ip yield in nucleus- nucleus 
collision as the function of the energy density or 
number of participants in the collision. The J/ip 
yield could show a significant suppression even if 
the energy density of the system is not enough to 
melt directly produced J/tp but is sufficient to melt 
the higher resonance states because they are loosely 
bound compared to the ground state J/ip. 

In nucleus-nucleus collisions, it is known that only 
about 60% of the observed J/tp originate directly in 
hard collisions while 30% of them come from the de- 
cay of Xc and 10% from the decay of tp'. Hence, the 

-integrated inclusive survival probability of J/tp 
in the QGP becomes [H,!!!]. 



iS" 



0.6(5' 



dir. 



0.3(5 



dir, 



0.1(5" 



dir. 



i-' 

(36) 

The hierarchy of dissociation temperatures in lattice 
correlator studies [53| (Table IV) thus leads to se- 
quential suppression pattern with an early suppres- 
sion of tp' and Xc decay products and much later 
one for the direct JM production. However, with 
our recent results [2J, [2^ (Table II-III) employing 
medium modification to the full Cornell potential 
and also results from potential model studies [56l | 
on dissociation temperatures, all three species will 
show essentially almost the same suppression pat- 
tern, i.e., the concept of sequential melting will not 
have any dramatic effect which will be seen in the 
next section. 



IV. RESULTS AND DISCUSSIONS 

There are three time-scales involved in the screen- 
ing scenario of J/tp suppression in an expanding 
plasma. First one is the screening time, Ts as the 
time available for the hot and dense system dur- 
ing which J/tp's are suppressed. Second one is 
the formation time of J/tp in the plasma frame 
{tp = jTp) which depends on the transverse mo- 
mentum by which the cc pairs was produced. Third 
one is the cooling rate which depends on the speed of 
sound, c^ through the equation of state. The screen- 
ing time not only depends upon the screening energy 
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TABLE II: Formation time (fm), dissociation tempera- 
ture Td (in units of 7^=197 MeV for a 3-flavor QGP) 
with the Debye mass in the leading-order [2^ . [25j , the 
speed of sound Cg and the screening energy density 
(GeV/fm^) calculated in SIQGP and ideal EoS for J/ip, 
ijj' , Xc states [2^.[25j. respectively. 



State 




Td 


(SIQGP) 


c?(Id) 


6, (SIQGP) 


e.(Id) 




0.89 


1.61 


0.26 


1/3 


17.65 


21.77 


^' 


1.50 


1.16 


0.24 


1/3 


04.51 


06.53 




2.00 


1.25 


0.24 


1/3 


06.15 


08.47 



TABLE III: Same as Table II but the dissociation tem- 
perature is calculated with the Debye mass in the lattice 
parametrized form [2^, [2^ . 



State 


TF 


Td 


(SIQGP) 


c?(Id) 


6, (SIQGP) 


e.(Id) 


J/^ 


0.89 


1.18 


0.24 


1/3 


4.83 


7.81 


^' 


1.50 


0.85 


0.18 


1/3 


1.21 


2.89 


Xc 


2.00 


0.90 


0.19 


1/3 


1.54 


3.36 



density, Cg but also on the speed of sound through 
EoS. The value of eg is different for different char- 
monium states and is calculated from the EoS under 
consideration and hence it varies from one EoS to 
other. If Eg > Ci, initial energy density, then there 
will be no suppression at all i.e., survival probabil- 
ity, S{pr^) is equal to 1. If the dissociation tempera- 
ture Td is higher, (oc T^) will also be higher, so 
system having > es will attain eg quickly. There- 
fore, the system will get less time to kill J/tp in the 
deadly region marked by the screening radius, re- 
sults in less suppression in the Jf-ip yield. However, 
for smaller values of T^, the system will take more 
time to reach resulting more suppression. More 
precisely, the screening time depends upon (i) the 
difference between the initial energy density and 
the screening energy density e^: the more will be the 
difference the more will be the suppression, (ii) the 
speed of sound: the values of which are less than 
1/3, the rate of cooling will be slower which, in turn, 
makes the screening time large for a fixed difference 
in (ei- Eg) leading to more suppression, and (iii) the 



TABLE IV: Same as Table II but with different values of 
dissociation temperature(s) obtained from lattice stud- 
ies O in units of rc = 175 MeV. 



State 


TF 


Td 


(SIQGP) 


c?(Id) 


6, (SIQGP) 


e.(Id) 




0.89 


2.10 


0.27 


1/3 


32.85 


60.83 


^' 


1.50 


1.12 


0.21 


1/3 


02.36 


05.80 




2.00 


1.16 


0.22 


1/3 


02.74 


06.53 



TABLE V: Dissociation temperature(s) Td (in units of 
Tc=203 MeV for a 2-flavor QGP) with the Debye mass 
in the leading-order [2^ . [25j of strong coupling. 



State 


TF 


Td 


(SIQGP) 


c?(Id) 


6, (SIQGP) 


es(Id) 




0.89 


1.69 


0.238 


1/3 


17.76 


23.14 


V.' 


1.50 


1.21 


0.172 


1/3 


04.07 


06.70 


Xc 


2.00 


1.31 


0.192 


1/3 


05.77 


08.79 



rj/s ratio: if the ratio is larger then the cooling will 
be slower, so the system will take longer to reach 
es resulting the higher value of screening time and 
hence more suppression compared to rj/s = 0. With 
this physical understanding we analyze (iS'(py)) as a 
function of the number of participants N-p^^^^ in an 
expanding QGP. 

In our analysis, we have employed the dissociation 
temperatures for J/ip, Xc a-nd ^p' recently computed 
by us m, m in the Table(s) II-III for 3-flavor and 
V for 2 flavor QGP, respectively, where wc employed 
an appropriate form of medium modified potential 
which includes both non-zero string tension effects 
beyond Tc and the usual screened coulomb term. To 
compare our predictions with the lattice correlator 
studies, we employ the dissociation temperatures (in 
Table IV) computed from the first principle of (lat- 
tice) QCD [5§| where a simple screened Coulomb 
form was assumed for the medium-dependent inter- 
quark potential. The corresponding values for and 
Cg calculated in the strongly-interacting and ideal 
EoS are also listed in the respective Tables. We shall 
employ these next to study {S(jj^)). 

We have shown the variation of -integrated sur- 
vival probability (in the range allowed by invariant 
spectrum of J/ip by the Phenix experiment [l|) 
with Np^j.^ at mid-rapidity in Figs. 1-8. The exper- 
imental data (the nuclear- modification factor Raa) 
arc shown by the squares with error bars whereas 

circles and diamonds represent with {{S )) with- 
out {{S )) sequential melting. We have employed 
three values of rj/s viz. 0, 0.08, and 0.3 in each fig- 
ure to the see the effects of dissipative terms in the 
expansion of the plasma. In Fig. 1, we find that 
the survival probability of directly produced J/ip 

is slightly higher than (S*^^*^^) and is closer to the 
the experimental results For the lower value of 
rj/s our predictions are closer to the experimental 
ones. As the ratio rj/s is increased, the expansion 
of the system becomes slower leading to the higher 
values for the screening time and hence lower values 
of (S'(py)) i.e., J/ip^s will be suppressed more. How- 
ever, for all the three values of 77/s ratio, {S{pj,)) for 
both the directly produced and sequential J/ip sit 
below the experimental numbers for Np^^^ 100. 
In Fig. 2, we used the same values of dissociation 
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FIG. 1: The variation of p^integrated survival probability (in the range allowed by invariant Pj. spectrum of J/tp 
by the Phenix experiment [l|) versus number of participants at mid-rapidity. The experimental data (the nuclear- 
modification factor Raa) are shown by the squares with error bars whereas circles and diamonds represent with 
incl dir 

{{S )) without {{S )) sequential melting using the values of Td's [241 ] and related parameters from Table II 
using SIQGP equation of state. 
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FIG. 2: Same as Fig. 1 but the related parameters from the ideal EoS. 
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FIG. 3: Same as Figure 1 but using using the values of Td's 
SIQGP EoS. 



and related parameters from Table III for the 
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FIG. 4: Same as Fig. 3 but employing the ideal EoS 



temperatures as in Table II but the thermodynamic 
variables viz. Cs, etc. have been calculated in 
the ideal EoS to see the sensitivity of the EoS to 
the plasma dynamics which, in turn, indirectly affect 
suppression in relativistic collision. In this case too, 
the 77/s dependence of {S{pj.)) pattern remains the 
same and (S'(py)) is always higher than the corre- 
sponding values (with SIQGP) in Fig.l. The match- 
ing is almost perfect for ri/s=0. This can be under- 
stood in terms of the sensitivity of {S{pj,)) on the 
speed of sound, for a fixed difference in (ei-Cs) be- 



cause cooling of the system with ideal EoS is much 
faster compared to the strongly-interacting EoS so 
the system will spend less time in the screening re- 
gion results in less suppression. 

Another interesting observation, which is common 
to both Fig.l and 2 is that as the ratio rj/s is in- 
creased from to 0.3, {S{p^)) for both directly and 
sequential J/t/i's become smaller. To be more ac- 
curate, (S'(py)) for the sequential decays will be af- 
fected much causing the suppression stronger and is 
being separated from the {S{pj,)) of the directly pro- 
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FIG. 5: Same as Fig. 3 but using the values of Tb's [53| and related parameters from Table IV for the SIQGP EoS. 



duccd J li). As the ratio rj/ s is increased from zero, 
cooling becomes slower so that Xc and ip' show more 
suppression due to their smaller value of tg (smaller 
Td) compared to J /ip making the difference between 
Ci and Eg larger. This leads to more suppression for 
X and V'"s. 

In Figs. 3 and 4, we used the dissociation tem- 
peratures employing the lattice parametrized form 
of the Debye mass (in Table III) which is the same 
as the values obtained from potential models but 
the expansion of the system is taken into account 
through SIQGP and ideal EoS, respectively. We 
again vary rj/s between and 0.3. Interestingly, 
in both the cases (Fig. 3 and Fig. 4) the suppression 
is large compared to the previous cases (Fig.l and 
Fig. 2) and is much smaller than than the experimen- 
tal results. Since the dissociation temperatures are 
much smaller than the values in Table II so that the 
difference (e^ — Cs) becomes large. Therefore, the sys- 
tem will take longer time to reach e^, as a result the 
system have enough time to kill J/ijj^s. The sensi- 
tivity to the EoS and the ratio 77/s is the same as in 
Figs.l and 2. As expected, there is more suppression 
in SIQGP as compared to ideal EoS. 

In Figs. 5 and 6, we use the dissociation tempera- 
tures from the lattice correlator studies [5^ (in Table 
IV) and the expansion of the system are taken into 
account through SIQGP and ideal EoS, respectively. 
Interestingly, in both Figs, the suppression is much 
smaller compared to the previous cases. The hier- 
archy in the dissociation temperatures in Table IV 
led the sequential suppression to play an important 
role which causes more suppression for the sequen- 
tially produced J/tp^s which is closer to the Phenix 
results compared to very less suppression for the di- 
rectly produced J/ip's. On the contrary, in the ear- 



lier sets on the dissociation temperatures the con- 
cept of sequential melting has not any dramatic ef- 
fect on (S'(py)). In this set of dissociation tempera- 
tures, the agreement between theory and experiment 
is better with the SIQGP EoS (Fig. 5) compared to 
ideal EoS (Fig. 6) because cooling is slower in the 
former EoS and results in more suppression for the 
sequential melting. In addition to the effect of EoS, 
the effect of the dissipative forces in terms of 77/s 
becomes prominent because the matching (with the 
experiment) is good for higher values of 77/s. This is 
due to the fact that the larger value of viscous force 
makes the expansion slower which results in more 
suppression. 

To examine the flavor dependence we have plot- 
ted the -integrated survival probability for the 2- 
flavor QGP with the dissociation temperatures in 
Table V in Figs. 7, 8. The suppression pattern re- 
mains the same in comparison to Figs.l and 2. It 
seems surprising because of the following reason: As 
the number of flavors in the system decrease, de- 
crease rapidly so that the difference (ci — es) becomes 
larger compared to 3-flavor system. As a result the 
system will get more time to kill J/i/j's resulting 
more suppression. But it did not happen in the Figs. 
7 and 8 because we take into account the flavor de- 
pendence of the dissociation temperature where the 
dissociation temperature for the 2-flavor system is 
much higher than 3-flavor sysem which compensates 
the decrease in the screening energy density due to 
the decrease in the degrees of freedom. This bal- 
ance makes the suppression pattern same in Fig.l 
and Fig.7 (or Fig. 2 and Fig.8). 

Let us now compare Fig.l with Fig. 5 and Fig. 2 
with Fig. 6. It is clear from Figs. 5 and 6 that {S{p^)) 
for both the directly and the sequentially produced 
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FIG. 6: Same as Fig. 5 but employing the Ideal EoS. 



J/ip are quite high with the higher values of Tp's 
obtained from the lattice correlator studies (in Ta- 
ble IV) compared to the values obtained from our 
model (in Table II). Since the dissociation temper- 
ature for the Jf-ip is high so that is much larger 
than Ci up to the number of participants 200 result- 
ing no suppression at all. However, T^i's are not so 
high for Xc and ip' so that the sequential probability 

{S""" ) has been suppressed that makes them closer 
to the experimental results for the higher values of 

rj/s. However, {S'"" ) gets 60% contribution from 
the directly produced J/tp, it has a dominating con- 
tribution on the behavior of sequential S'(p^). This 
is indeed reflected in Fig. 5 and Fig. 6. To be conclu- 
sive, the dissociation temperatures exploited in the 
present work (in Table II) shows better agreement 
with the experimental results from RHIC as com- 
pared to the spectral function technique calculated 
either in the potential model [5^ or in the lattice 
(temporal) correlator [58l |. 



V. CONCLUSIONS AND FUTURE SCOPE 

In conclusion, we have studied the charmonium 
suppression in a longitudinally expanding QGP in 
the presence of dissipativc forces. Wc find that pres- 
ence of dissipative terms in the fluid equation of mo- 
tion slower the expansion rate and eventually lead 
to the enhanced suppression of J/ip. In other words, 
the presence of viscosity enhances the screening time 
for J/tp in the SIQGP medium and hence the sur- 
vival probability gets decreased compared to that 
without the viscous forces. These conclusions are 



true for both the directly and sequentially produced 

J/iK 

In this work, we have exploited a recent under- 
standing of dissociation of quarkonia in the QGP 
medium which rely on the fact that the transition 
from the hadronic matter to QGP is a crossover not 
a phase transition in the true sense. We have em- 
ployed the results of [23, on dissociation tem- 
peratures of various charmonium states. We have 
employed the SIQGP equation of state to estimate 
the screening energy density and the speed of sound 
to study the J/ijj yield. To compare our results with 
those obtained by employing the simple screening 
picture of quarkonia commonly considered in the 
literature, we employ dissociation temperatures esti- 
mated in lattice correlator studies [sHi ■ We find that 
the results on J/ip survival probability agree with 
the Phenix Au-Au data[l| with the set of dissocia- 
tion temperatures (Table II) obtained with the per- 
turbative result of the Debye mass. We had shown 
in our earlier work [2^ that the inclusion of non- 
perturbative contributions to the Debye mass lower 
the dissociation temperatures substantially which 
looks unfeasible to compare to the spectral anal- 
ysis of lattice temporal correlator of mesonic cur- 
rent [5^ which finally makes the J/^ survival prob- 
ability too small to compare with the experimen- 
tal results. This does not immediately imply that 
the non-perturbative effects should be ignored. It 
is rather interesting to investigate the disagreement 
between the non perturbative result obtained with a 
dimensional-reduction strategy and the Debye mass 
arising from the Polyakov-loop correlators. Only fu- 
ture investigations may throw more light on this is- 
sue. 
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FIG. 7: Same as Fig.l but using the values of Td's and related parameters in SIQGP EoS from Table V for 
2-flavor QGP. 



This leaves an open problem of the agreement be- 
tween these two kind of approaches. This could 
be partially due to the arbitrariness in the crite- 
ria/definition of the dissociation temperature. To 
examine this point we estimated both the upper and 
lower bound on the dissociation temperatures [13, 
[25l |. Thus, this study provides us a handle to deci- 
pher the extent up to which non-perturbative effects 
should be incorporated into the Debye mass. 

Finally, our study does provide a systematic way 
to analyze the fate of charmonium in QGP medium. 
We have included two important aspects of the QGP 
at RHIC which are based on recent experiment find- 
ings: (i) strongly- interacting picture of QGP and (ii) 
non-vanishing string tension (between qq) contribu- 
tions beyond the deconfimcnment point. The first 



one we have incorporated by considering the phe- 
nomenological equation of state [3^ and the speed 
of sound determined from it and the second one is 
by considering a right criterion for the dissociation 
of quarkonia in QGP. Our attempt is perhaps the 
first one to understand J/ip suppression systemati- 
cally in SIQGP. It would be of interest to extend the 
present study by incorporating the higher order con- 
tributions coming from the viscous forces including 
contributions of the bulk viscosity. These issues will 
be taken up separately in the near future. 
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